Large Cross-Phase Modulation between Slow Co-propagating Weak Pulses in "^^Rb 
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We propose a scheme to generate double electromagnetically induced transparency (DEIT) and 
optimal cross phase modulation (XPM) for two slow, co-propagating pulses with matched group 
velocities in a single species of atom, namely *^Rb. A single pump laser is employed and a homo- 
geneous magnetic field is utilized to avoid cancelation effects through the nonlinear Zeeman effect. 
We suggest a feasible preparational procedure for the atomic initial state to achieve matched group 
velocities for both signal fields. 

PACS numbers: 42.65.-k, 42.50.Gy, 03.67.-a 



Introduction: — The optical Kerr effect, n = n^^^ + 
rS^^I (for n the total refractive index, n'^") the linear re- 
fractive index, / the field intensity, and n*^^-* the optical 
Kerr coefficient) is invaluable for spectral broadening and 
self- focusing of laser pulses . Large nonlinear interac- 
tions have also been used to generate single photons 3, 
enhanced refractive index 0], self-phase modulation ^, 
and Stokes and anti-Stokes generation 0, . In cross- 
phase modulation (XPM) , / is the intensity of the other 
field; thus XPM enables the phase of each field to be 
controlled by the strength of the other, which is criti- 
cal for applications such as deterministic optical quan- 
tum computation and all-optical switching [^. Un- 
fortunately vP'^ is extremely small but can be effectively 
increased via electromagnetically induced transparency 
(EI T) Ip. Iiol llll . w hich yields a large optical nonlinear- 
ity list Il4 Il5l|. Whereas these experiments demon- 
strated the promise of a large nonlinear coupling for cw 
fields, applicability is severely limited for propagating 
fields due to the different group velocities of the two 
pulses, which significantly reduces the interaction time. 

Double EIT (DEIT) generalizes EIT for simultaneous 
action on two separate fields, and, for slow group veloc- 
ities, can effect a large XPM HI [13, III Ull; however 
DEIT has not yet been successfully realized. In this let- 
ter we devise a method for achieving large XPM for two 
slow co-propagating pulses with matched group veloci- 
ties in a single species of atom, namely using the Dl line 
of ^^Rb, by avoiding cancelation of nonlinearities near 
resonance and decoherence due to coupling bright states, 
plus employing population transfer to achieve matched 
group velocities for the two pulses. 

Scheme: — Our scheme employs five atomic levels as 
shown in Fig. QJa) and combines the advantageous prop- 
erties of the N-type scheme the tripod scheme 0| 
and the M-type scheme A single pump field drives 
the 1 3) ^ 1 4) transition with Rabi frequency fip. If 
state 1 2) is omitted our scheme reduces to the N-typc 
scheme: the AC Stark shift created by off-resonant cou- 
pling to state 1 5) creates a giant nonlinearity between 
signal fields 1 and 2, but EIT is only realized for signal 
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FIG. 1: Atomic level configuration for creating DEIT and 
large XPM between two weak signal fields 1 (fii) and 2 
{yi2,0!2)- (a) Simplified level configuration of the scheme, 
(b) Realization of the scheme using the Dl Line of ^'^Rb. 



field 1. Consequently, due to the mismatch in the group 
velocities, the interaction time is greatly reduced. On the 
other hand, if state |5) is omitted we have a tripod scheme 
in which DEIT can be achieved but for which the non- 
linearity disappears at exact two-photon resonance (as it 
does in the M-type scheme). The slight detuning that 
is therefore necessary to create nonlinear effects 0, 0| 
also generates linear absorption, and for a given desired 
absorption rate the generated nonlinearity is suppressed 
by one to two orders of magnitude. 

We will show that the scheme of Fig. ^a) avoids the 
problems associated with the schemes discussed above. 
In addition to the pump field, two signal fields with slowly 
varying field amplitude £i {i ~ \,2) resonantly couple the 
|1) ^ |4) and |2) |4) transitions with Rabi frequencies 
51i = — |d4i|£i/?i, with dy = (i|d|j) the matrix elements 
of the dipole moment operator d. In addition, field 2 
also off- resonantly couples the |3) ^ |5) with Rabi fre- 
quency ^,'2 = — |dR.ci|£^//l. If the atoms are initially pre- 
pared in a mixture [20| of state |1) and |2) with density 
matrix pniix = +P2|2)(2|, the two A-subsystems 

|1) ^ |4) ^ |3) and |2) ^ |4) |3) induce EIT for 
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both fields 1 and 2, i.e., DEIT. The real part of the cor- 
responding index of refraction Ui for signal field i then 
has the approximate form [loj 

n, = 1 + r]iSp^ , ?7, = ppi |d4i| V(2feo|f^pn (1) 



for i = 1,2, 5ij = Sj — Si, p the atomic number density, 
and Pi the initial population in state For Eij = Ej — 
Ei, the detuning of the signal field frequency uJi and the 
laser pump frequency ujp are Si = tOi ~ Ei^^/h and 5p = 
cijp — E^i/h, respectively. Because Eq. ^ varies strongly 
with the two-photon detuning Spi the group velocity is 
very small. Below we will show how the group velocities 
can be made equal by manipulating the populations pi. 

The off-resonant coupling for transition |5) ^ |3) 
caused by field 2 modifies the tripod scheme by producing 
an AC Stark shift for state |3) given by 



AE3 = h\n'^\'/A, A = uj2-E35/h 



(2) 



with A the (large) detuning of field 2 with respect to 
the |3) <-> |5) transition. This energy shift implies that 
the pump frequency ujp is somewhat detuned from the 
1 3) <-> 1 4) transition; hence the index of refraction for 
this EIT medium changes due to the dispersion relation 
of Eq. (HJ. 

Replacing Spi by Spi — Ai^a/fi, and inserting this into 
Eq. iQJ yields 

ni = 1 + 771 (5pi - xh), X = |d53/?i|Vc£oA (3) 

with intensity /i = ^E^xB*! = ceolifiP corresponding 
to the modulus of the complex Poynting vector. 

We find that this intuitively appealing phenomenolog- 
ical derivation agrees with a rigorous calculation based 
on third-order time-dependent perturbation theory in the 
weak signal fields. Neglecting all decoherence effects and 
assuming a resonant pump field {Sp = 0), the Schrodinger 
equation yields the index of refraction 



"-1 = 1 + Vii^i - Xh) 
"-2 = 1 + ri2{S2 - xh) 



(4) 
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with the self-phase modulation (SPM) coefficient for 
field 2 given by r]2X and cross-phase modulation (XPM) 
coefficient for both fields given by rjix- Both XPM and 
SPM coefficients are given to first order in 1/A and to 
second order in l/|r2pp. A more detailed analysis (see 
below) shows that this simple result is qualitatively cor- 
rect if the detuning A is much larger than the decay rate, 
and the AC Stark shift (O is smaller than the width of 
the EIT transparency window. As this scheme shares the 
advantages and avoids the disadvantages of the N-type, 
tripod and M-type schemes, it yields the optimal XPM 
based on EIT techniques. 

Implementation with Rb: — Our scheme can be re- 
alized in a gas with a single species of atoms; a specific 



implementation using the Dl line of *^Rb is illustrated in 
Fig. nfb) . A homogeneous magnetic field parallel to the 
laser propagation minimizes coupling to states that are 
not part of the scheme in Fig.^^a). It is of particular im- 
portance to break two-photon resonance for the A subsys- 
tem |3) ^ |5) ^ |X). If this is not the case a further EIT 
transition that transfers atoms from state |3) to state |X); 
this would remove the crucial energy shift ^ and hence 
destroy XPM. The linear Zeeman effect is not suitable to 
break two-photon resonance. However, for large enough 
magnetic fields the Zeeman splitting depends nonlinearly 
on the magnetic quantum number without affecting the 
corresponding selection rule |2ll | . Numerical diagonaliza- 
tion of the Hamiltonian for a magnetic field of i? = 150 G 
shows that the energy differences E32 and Exs differ by 
an amoimt (5mag = {E32 - Ex3)/h = -12.9 MHz. This 
is much larger than the EIT transparency window 
St ~ min{r2p, |f^pP/7} and breaks two-photon resonance 
for the 1 3) <-> |5) <-> |X) transition. We therefore can ne- 
glect state |X) in our considerations; this approximation 
is confirmed by our numerical simulations (see below). 

We analyzed this specific scheme in three ways: the 
simple analytical theory (SAT) described above yields 
Eq. for the XPM coefficient, a more elaborate ana- 
lytical theory (EAT), and a numerical simulation (NUM). 
The elaborate theory EAT is based on third-order time- 
dependent perturbation theory to solve the Schrodinger 
equation with the evolution restricted to states |1) ^ |7) 
of Fig. n^b). Spontaneous emission is included in the 
non-Hermitean atomic Hamiltonian 
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with detunings (see Fig. ^b)) Ai ~ uJi — Eiq/H, A2 = 
UJ2 — E^t/H and phenomcnological complex quantities 



Ai ~ — Ai — i^ [23. We obtained an accurate analytical 
result, but the expression is unwieldy and not displayed 
here. Instead we compare it with the numerical simula- 
tion NUM in which we solved the time evolution of the 
atomic density matrix based on the full Lindblad master 
equation for all 16 states of Fig.njb), including sponta- 
neous emission based on all Clebsch-Gordan coefficients 
and 6-j symbols as specified by Steck Below we 

will also briefly discuss numerical results that include de- 
phasing and transit broadening which is important for 
hot atomic gases. XPM between the two signal beams is 
proportional to that part of the atomic mean dipole mo- 
ment which couples to signal field i and vanishes if the 
intensity of the other signal field is zero. To describe this 
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we define a dimensionless mean cross dipole moment by 
dxPM = ^ (|Tr(pd)|^^__^o - |Tr(pd)|^^_^^o) , (6) 

with e the electron's charge and ao the Bohr radius. 
dxPM is a direct measure of the mutual inteaction be- 
tween the signal pulses; in SAT and EAT it is propor- 
tional to the XPM terms, in NUM it also contains higher 
order terms. The result of NUM for dxPM are shown in 
(Fig.|2J) and show excellent agreement with the analytical 
results EAT and qualitative agreement with the simple 
theory SAT. We have used Steck's spectroscopic data |23 
and the parameters given in the figure caption. 

The XPM phase shift can be identified as the difference 
between the phase factors acquired by signal field i if the 
second signal field is switched on and off, respectively, 
after it has been propagating for L = 1.6 mm (double 
the Rayleigh length). This definition amounts to 

^XPM ^ ^^l(„,^(£3_^ ^ 0) - n,(53-. = 0)) , (7) 

for i — \,2 and with the refractive index evaluated us- 
ing the result of EAT for the steady state value of the 
atomic dipole moment. Fig. |5] indicates that steady 
state is already achieved after 0.7 /is. The stationary 
values are given by (fi^^^ = 1.11, corresponding to an 
effective XPM coefficient {-qixf^ = 0.27cmVW, and 
^XPM _ Q_g4^ corresponding to an effective XPM co- 
efficient (r?ix)2^ = 0.20cmVW. The difference of XPM 
for field 1 and 2 is due to dccoherence effects and cou- 
pling to off-resonant states. The pulses are attenuated 
by about 20% due to off- resonant coupling to states |6) 
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FIG. 2: Mean cross dipole moment dxPM of Eq. Q as a 
function of time for signal field 1 and 2. Lines 1, 2 represent 
the result of SAT (Eq. Q), which is identical for both signal 
fields. Lines 3, 4 (5, 6) represent the results of EAT and 
NUM for signal field 1 (2), respectively. The parameters used 
for these solutions are Q.i = 0.68 MHz, ^2 = -0.55 MHz, 
Q.p = 4.06 MHz, A = -134.58 MHz, Ai = 894.93 MHz, 
A2 = 621.85 MHz, 5mag = -12.93 MHz, and 7 = 5.73 MHz, 
p = 10^* cm~"^ and B = 150 G, with initial atomic state 
Pmix = 0.4|l)(l|-h0.6|2)(2|. 



and 1 7). In addition, the nonzero detuning associated 
with the energy shift (01 leads to an intensity-dependent 
attenuation of about 3%. 

In summary, the three theoretical results SAT, EAT, 
and NUM indicate that the five-level scheme of Fig.^b) 
would work well even if all 16 hyperfine states arc in- 
cluded. The off-resonant coupling to the additional states 
only results in the small, very rapidly oscillating pertur- 
bations displayed in Fig. [5] as an increased width of lines 
4 and 6. 

Preparation of the initial state: — For pulsed signal 
fields maximum interaction time is achieved for DEIT 
with equal group velocities. From Eq. I^J one can in- 
fer that the group velocity Vi of signal field i is propor- 
tional to Vi oc Pi|d4ip. It is therefore possible to achieve 
f 1 = V2 by preparing suitable populations pi of the ini- 
tial density matrix pmix = pi|l)(l| -|-p2|2)(2|. We pro- 
pose the following procedure. We first pump all atoms 
to state |4) whence they decay to states |1) and |2) plus 
other states. To maximize the EIT effect one may re- 
pump atoms in all states but |1) and |2) into state |4) 
until the atomic state is well approximated by Pmix with 
Pi = 74i/(74i + 742) oc |d4ip, where 74, is the decay rate 
from state |4) to state With these populations the 
group velocities would be proportional to Vi cx |d4i|'* and 
thus differ significantly. We therefore induce as a last step 
in the atomic state preparation a Raman transition be- 
tween |1) and |2). Ideally this exchanges the populations. 
Pi <^ P2, resuhing in p, = |d4,3_ip/(|d4ip| -h d42p) 
so that both group velocities are now proportional to 
|d4ip|d42p and therefore equal. We remark that the 
Raman transition is not required to work perfectly. An 
exchange of only 90% of the populations would lead to a 
difference of 10% in the group velocities, for instance. 

Maximum phase shift for pulses at single-photon 
level: — One prominent application of XPM would be the 
creation of a controlled phase gate for photonic qubits 
[2^. We therefore estimate the maximal XPM phase 
shift achievable for two Gaussian pulses at single-photon 
level propagating through DEIT media. To do so we gen- 
eralize the expression of XPM Kerr coefficient by re- 
placing 7/1 by the standard EIT expression that includes 
spontaneous emission |^ for A atoms in a resonant 
pump beam, t/i = ?7i|r2p| V(|^^pP - h{h + i"f/2)). To 
first order in the detuning Si we have fji w 771 -I- irSi with 
T = 7/i7/(2|r2pp). Within the paraxial approximation, 
the slowly varying electric field amplitude of a Gaussian 
pulse of duration T and minimum (1/e intensity) waist 
wq at single-photon level in a medium that is described 
by fji propagates according to 

s = ^---"es — > ^ ^ 

/ ^ 
V V2^eocTwl 
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with Q : 



and ^ = + 4fczr. As for any Gaus- 



sian beam the field intensity is only close to its maximal 
value over a propagation length L = 2zr, with = kwQ 
the Rayleigh length. In addition, the EIT transparency 
window attenuates the pulse through the factor S. Re- 
quiring that S changes little over zr implies a minimum 
2\/Tkzii. Using these values we 
^ as the prod- 



pulse duration of T 
can estimate the maximum phase shift 
uct of L and the XPM induced change in the wavevector 

k = nuj/c, 



L-x 



(10) 
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To avoid decoherence we need 7/ A <C 1. The diffraction 
limit |l7| implies X/wq < 1 and we estimate that res- 
onant dipole-dipole interaction would modify (j)"^^^ for 
pk^^ ~ 1. Thus, the maximum XPM phase shift be- 
tween two photons would in this approach be of the or- 
der of 0.1 rad. This is considerably smaller than previ- 
ous estimates of more than one rad 0, 0| but may be 
overcome by combining our scheme with the ideas pre- 
sented in Ref . . The difference arrives from imposing 
different constraints (dipole-dipole interaction) and our 
neglect of ground-state decoherence, which is reasonable 
for trapped ultracold atoms. 

Limitations: — The numerical example presented 
above is based on the parameters for a dense ultracold 
gas such as a trapped elongated Bose-Einstein conden- 
sate. A mutual phase shift of the order of tt could then 
be obtained for pulses containing hundreds of photons. 
In a hot atomic gas (e.g., T=440K and p ~ lO^^'cm"^), 
the efficiency of this scheme is mainly limited by the 
atomic motion. Although the Doppler effect is canceled 
in the two-photon detuning 6i — 6p for co-propagating 
light fields, the atomic motion narrows the transparency 
window according to St < ^Ip/Ao where A/^ « 

500 MHz is the mean Doppler detuning. The XPM phase 
shift would then be two orders of magnitude smaller than 
the optimal value. Transit broadening and ground state 
dephasing due to the atomic motion in and out of the 
laser fields and collisions between atoms are also detri- 
mental. For laser pulses with a width of 2mm in a hot 
atom setup using buffer gas 0], rates less than IkHz are 
possible. Numerical simulations indicate that our scheme 
would produce 0.4rad XPM (with 40% attenuation) for 
two weak pulses with 10^ photons, a transit broadening 
rate of IkHz, and a ground state dephasing rate of IkHz. 

Conclusion: — Cross-phase modulation between weak 
signal fields is one of the most significant challenges for 
nonlinear optical switching and quantum information. 
We have proposed an approach to achieve strong cross- 
phase modulation under realistic conditions for *^Rb in 
which matched group velocities for two interacting pulses 
can be realized. Our approach incorporates all of the fol- 



lowing desirable conditions: co-propagating laser beams 
to avoid Doppler shifts, a single atomic species, and only 
a single pump field. We show that our scheme yields op- 
timal XPM for pump schemes based on DEIT, and we 
estimate the maximum phase shift to be 0.1 rad for Gaus- 
sian pulses at single-photon level. To avoid decoherence 
effects Bose-condensed gases would be preferred, but we 
expect the scheme also to be effective for thermal gases 
if sufficiently strong signal pulses are used. 
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